R=120;
L=100e-6;
tau=L/R;
T=tau/100;
t=0:T:(100*tau);
N=length(t);
I(1)=0;
E=1;
for k=2:N;
    I(k)=I(k-1)-(3*T*R/(L))*I(k-1)+(E*T)/(L);
end;
plot(t,I*R);
